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We consider Einstein-Maxwell-Dilaton (EMD) Lagrangian supplemented by double Liouville po- 
tentials to enrich our system and investigate the resulting dynamics. The general solution provides us 
alternative structures induced on the 3-dimensional domain wall (DW) moving in the 4-dimensional 
bulk. In particular, the local maximum in the potential suggests a maximum bounce (or onset for a 
contraction phase) of the 3-dimensional Friedmann-Robertson- Walker (FRW) universe on the DW. 
Depending on the choice of parameters we plot various cases of physical interest. 
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The idea that our universe is a brane living in a higher dimensional space has attracted considerable attention during 
the recent decade |f|. The physical properties on such a brane must automatically be induced from the surrounding 
space known as bulk through some well - established junction rules. One such particular example of brane - bulk pair 
is provided when the brane has dimension one - less (d — 1) in a d— dimensional bulk known as domain wall (DW) 
with the Z2 symmetry across. Beside the universe concept, the DW entails topological defect formations pertaining 
1— — to the remnants from the big bang, different vacua of quantum fields and other applications of vital importance. 

As emphasized often, one of the most important problems in modern cosmology is the accelerated expansion of our 
^ ' currently observed universe, its causes and whether it will do so indefinitely. Can the answer be obtained from the 
idea of DWs assumed that our universe is a DW in a higher dimensional bulk spacetime? In other words, can we 
construct an appropriate bulk spacetime so that inflation driven by matter on the DW is induced from it? It must 
be admitted that a bulk with such generality is still missing. Herein we choose our bulk metric to be 4— dimensional 
£\| \ (d = 4) while the DW becomes 3— dimensional (d — 3) which may be considered in this regard as a toy model. 
^ 1 Let us not forget, however, that lower dimensional physics may serve to shed light on higher dimensions. One 
justifiable reason for choosing our bulk as d = 4 is that we can provide an exact solution with quite generality to 
encompass new dynamics on the DW and serve our purpose well. The action contains Maxwell and dilaton fields 
beside gravity whose coupling is non-minimal 0, Q . With this much rich action, unfortunately, we were unable to 
obtain exact solution in bulk dimension d > 4, for this reason we are satisfied with d = 4. In addition, we have an extra 
Liouville type potential of the dilaton Q and our Maxwell field is pure magnetic. Due to the fact that the coupling 
to dilaton is rather complicated, solution becomes possible, with a pure magnetic field. Depending on the choice of 
integration constants and dilatonic parameter we have obtained a large class of non-asymptotically magnetic solutions 
that yield previously known solutions of its kind in particular limits • Pure electrically charged solutions Q , on 
the other hand are not expected to overlap with our magnetic ones. One crucial difference between those and the 
present work is in the coupling between dilaton and the magnetic field. Essentially it is this difference that provides a 
richer dynamical structure on the DW without reference to non-physical conditions such as negative mass. It should 
also be added that we consider a neutral brane (DW) so that by Z2 symmetry, continuity of the vector potential and 
electromagnetic field across the brane the Maxwell equations are trivially satisfied. In other words the surface (DW) 
Lagrangian doesn't depend on the electromagnetic field, it depends explicitly only on the dilaton whose boundary 
conditions are accounted for in the junction conditions. This amounts to the fact that the tension / energy density on 
the brane are only those induced from the energy and dilaton in the bulk. The effect of the Maxwell fields manifests 
itself on the brane through the magnetic charge. 

Our DW universe is a 3— dimensional Friedmann - Robertson - Walker (FRW) universe with a single metric function 
(=the radius of the universe) depending on its proper time Q. The boundary conditions that connect bulk to the 
DW are provided by the Darmoise - Israel conditions Q apt for the problem. These conditions determine the energy 
- momentum on the DW universe together with its dynamics. An interesting aspect of the solution that we present 
in this paper is that on the DW our FRW universe attains both lower and upper bounces. This particular point 
provides our main motivation for this study. It happens that the occurrence of the double Liouville type potentials 
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in our model made this possible for us. From the boundary conditions the induced potential on the DW the radius 
a (t) of the DW universe satisfies an equation of a particle with zero total energy (i.e. (J 2 ) + U (a) — 0). Such a 
motion is physical only provided the potential U(a) satisfies U(a) < 0, for all a(r). This is investigated thoroughly 
and we observed that our DW universe admits both minimum and maximum bounces. Among other things absence 
of a DW universe (i.e. U > 0) is also a possibility. 

In Section II we present our model Lagrangian and obtain a general class of solutions to it. In Section III we analyze 
the induced potential on the DW. Conclusion makes our last section IV which summarizes our results. 



II. OUR MODEL IN d = 4 AND ITS SOLUTION . 

Our action of Einstein-Maxwell Dilaton (EMD) gravity is written as S 

K JM \Z Z A J k Jj2 Js 

where T = F Xa F Xa is the Maxwell invariant, V ((f)) = V x e p ^ + V 2 e ^, W ((f)) = Aie -2 ^ + A 2 e" 27 ^ and = (r) is 
the dilaton scalar potential. Herein j i , f3 i , Vj and \ are some constants to be identified later while Cdw = — V (4>) = 
—Voe e ^ (e — const.) is the induced potential on the DW. {K} is the trace of the extrinsic curvature tensor of 
DW with the induced metric hij (h = \gij\). (Latin indices run over the DW coordinates while Greek indices refer to 
the bulk's coordinates. Also in the sequel we use units in which k 2 = 8irG = 1). The 4— dimensional bulk metric is 
chosen to be 



1 

m 



ds 2 = -f (r) dt 2 + -yj^dr 2 + H(r) 2 dn 2 , (2) 



where / (r) and H(r) are functions to be found and dQ| is the line element on a 2— dimensional space of constant 
curvature with k € { — 1,0, +1} , i.e., 

C d0 2 + sm 2 9dip 2 , k=l 
dQ 2 k =l d6 2 + dip 2 , k = . (3) 
[ dd 2 +sinh 2 6d(p 2 , k = -l 

The field equations inside the bulk follow from the variational principles as 

R^=d^d^ + V(^)S; + W(4>)T», (4) 

v 2 = y'(0) + ^'(0)j-, (5) 



in which 

T; = 2F, X F ^ l -F&l (6) 

is the energy momentum tensor of our Maxwell 2-form F =\F VLV dx il A dx v . Variation with respect to the gauge 
potential 1— form A yields the Maxwell equation 



d (W (0)* F) = 0, (7) 

in which d(.) is the exterior derivative and the hodge star * means duality. As we commented before, is a function 
of r, and W ((f)) — Aie -271 ^ + A2e -272 * is double Liouville function. This means that although an electric field ansatz 
that makes the Maxwell equation complicated enough for an exact solution, a pure magnetic ansatz easily satisfies the 
Maxwell equation. Therefore we prefer to use a magnetic potential 1— form with magnetic charge P which reads as 
(This is also another feature of our work in comparison with Maity's work in Ref. [5j, while in the case of Yazadjiev 
p both electric and magnetic fields are used) 

!—Pcos8 dtp, k = 1 
P 9 dip, k = . (8) 
Pcoshfl dip, k = — 1 
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The field 2— form becomes 



Psin9 d8A dip, k = l 
P dO A dip, k = 
Psinhe dO A dip, k = -1 



(9) 



where the *F is substituted into Eq. (7) and it is easily satisfied which justifies at the same time our choice of a pure 
magnetic gauge potential. For the case of k = 1, the magnetic charge p is defined as 



P 



— i F P 

Air 

S 2 



(10) 



where S 2 is the two-dimensional sphere while for the cases k = 0, — 1 the topological charge takes the form p = 

in which to is the area of the corresponding 2— surface. Based on our choice of the magnetic field, one can show that 

independent of k the Maxwell invariant reads 



and the energy momentum tensor becomes 



2P 2 
H(r) 4 



1 



-T diag [-1,-1, 1,11 



(11) 



(12) 



Upon substituting these into the field equations (4-5) we get the following equations 



— [f'H + 2f'H'} 
i [f'H + 2 f'H' + AfH"} 
±[HH'f + fH' 2 + fHH"]-^ 



= -V + -WF, 

= -f$ 2 -V+ X -WT, 

= -V-\WT, 



(13) 
(14) 
(15) 
(16) 



We note that a prime " "' implies derivative with respect to r and d,p = ^ . Before we write the solution of the field 
equations we comment that our interest here is in finding an exact, non-asymptotically flat solution which partially 
was found previously [lfj. Whether similar attempts may give a different class of asymptotically flat solution will 
remain as an open problem. 

The general solution for the metric function (after setting /3 1 — a\/2, (3 2 — „ > 7i = — ^ an d 7 2 = ) is expressed 

by 



f(r) 



1 + a 2 



2(2 + c 2 ) 

A 4 (l+a 2 ) 



a 2 (l+a a ) 



2M (a) 



3 + o 
1 + n 



(tJ + 



2 ■ 4 / P 2 A 2 
3 I - 4^" 



-F 2 -^(ajhU -32M(^ 



a 2 -3^0 



a 2 -3 = 



(17) 



Other functions read as 



Jo. 



1 + a 



where the constant A is related to P by the constraint condition 



V 2 (l- a 2 ) A 4 + ka 2 A 2 - X 2 P 2 (1 + a 



0. 



(18) 



(19) 



We also note that in the solution (17) M and r are two integration constants. We shall use this general solution, 
with k = 1, in the following Section to construct our 2 + 1— dimensional DW. Solution (17) has five model parameters 



4 



(Ai, A 2 ,Vi,V2 and a) and three free parameters (after considering the constraint (19)) (P, M and r ) which gives us 
a large class of BH or non-BH solutions. To get closer to the solution we give some limits which may be useful in 
future calculation. First, we consider < a 2 < 1 and we find the asymptotic behavior of the master solution as: 

1„/ (l + a *)V(-^(^)*) =s M -Vi)oo, (20) 

lim/ (1W)V Lf^a) (7)^) = »P»M», (21) 

which show that by choosing proper values for V\ and Ai we definitely will obtain BH cases, at least with single 
horizon. Second, we consider a 2 = 1 which in the extremal limits admits 



lim / — > 2rro 



P 2 \ 



lim f -> 2r 

r-yO 



.4 4 

a 

2 



2 - V 2 - Vi 



= *H«[^-V2-Vi)>:. (22! 



=,^0 oo. (23) 



The third case is the choice of a 2 > 1 which has the limits as 

^ {-Ar- V2 )[T) ='!>»[-&- ( 24 ) 

lim/ = sgn (Ai) oo. (25) 
r— >0 

Here also in a similar manner by setting proper values for Ax, A 2 and V 2 we can construct BH with at least one horizon. 
(We note that the limit of the metric for r —> in all cases are the same.) Among interesting cases one may find that 
the choice - V 2 = gives different asymptotic limits which for < a 2 < 3 are given by 

2\2 2 I V l f r 0\~c 



lim / -^"T^^^J^ ) = '9n{-Vi) oo, (26) 

lim/ = sgn (Ai)oo, (27) 
r— >0 

while for 3 < a 2 

3+t» 2 

lim / -2M (l + a 2 ) 2 r 2 f — ) T7 " lr = s. 9 n (-2M) oo, (28) 
lim/ = sgn (Ai)oo. (29) 

In both cases we have the possibility of choosing proper values for the free parameters to have BH with at least one 
horizon. It should also be added that, to construct a DW in these bulk solutions we choose the radius of the DW 
always larger than the possible event horizon. This guarantees that the DW will not face any singularity problem on 
its domain. 

Let us add that, previously reported solutions are recovered when some constants of the present model are set to 
zero, e.g. A2 = V% = leads to the type II solution in Ai = V\ — leads to the type III solution in [4], etc. 

After all these asymptotic consideration we investigate the form of singularities at r = 0. To do so first we set ro = 1 
and obtain the Ricci Scalar as 

f 4A 2 (« 4 +a 2 + l)V2-fc( a 4 + l) 6(a 2 -2)V! 4a 2 M , 2P 2 a 2 A 1 ^2 / o 

R= ) A2 (1 + a 2 ) 2 r 2/(l + ^) + ( Q 2_ 3)l .2 Q 2/(l + o2) (1 + Q 2 ) 2 r (3 + Q 2 , / (l + Q 2 , + ( ^ " — ( 4+2 „2 ) y (l + 2 J , 

I 26y 2 ^ 2 ~5fc _ 6M-32Vi+12Vi lnr , 3AiP 2 2 _ o 

(30) 

It clearly shows that in any possible case the solution is singular at r = 0, which for BH cases are screened by 
horizon(s). 
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III. INDUCED POTENTIAL ON THE DW AND ITS IMPLICATIONS 

The 3— dimensional DW on the surface E in a 4— dimensional bulk A4 splits the background bulk into the two 
4— dimensional spacetimes which will be referred to as Ai±. Here ± is assumed with respect to the DW. Let us look 
at the master solution (17) and its asymptotic behaviors given by (20) - (29). We set the parameters such that 
lim r _ i . 00 / (r) = oo and for the case of BH we choose < r = a. For the non-BH case we make the choice < r = a. 
Upon imposing the constraint 

with the DW at r = a (r) , r being the proper time with respect to the wall observer, the DWs line element takes 
the form 

ds 2 dw = -dT 2 +a{T) 2 dn 2 k . (32) 

This is the standard FRW metric in 3— dimensions whose only degree of freedom is a(r), the cosmic scale factor. 
Now, we wish to employ the general solution for the field equations (13)-(16) under (11) and metric ansatz (2). We 
impose now the rules satisfied by the DW as the boundary of M.±. These boundary conditions are the Darmois-Israel 
conditions which correspond to the Einstein equations on the wall Q . These conditions on the DW E are given by 

-(J K Ki)-{K)Si)=Si, (33) 

where the surface energy-momentum tensor Sij is given by [3, d, H| 

Sv = ^^f#W=K(-n*j). (34) 
Note that a bracket (.) , in (33) implies a jump across E. The stress-energy tensor reduces to the form 

Si = -V(<f>)S>. (35) 

By employing these expressions through (33) and (34) we find the energy density and surface pressures for generic 
metric functions / (r) and H (r) with r = a (t). The results are given by 



(36) 



S$ = S*=p e =p v = ( /' + . 2a . 2 + 2y/f(a) + a*%) (37) 

V v / ( a ) + « H ) 

in which a dot "•" and prime ""' means -4- and 4-, respectively. The Einstein-equations on E accordingly read 

y/fJjp^E = Av(4>), (38) 

f i On ff' 

J - +2V/(a) + a 2 — = -V (</>). (39) 



We observe first that 



which after using (38) it becomes 



f' + 2a 2d 



y/f (a) + a 2 ddr 



(v//(«) + « 2 ) (40) 



y/f (a) + a 2 ddr\ 4 ^ ' H' J da \ 2 yT ' H 



6 



Finally as a result with (39) and (38) yields 



V{4>). 



da V H' 

This equation admits a simple relation between H(r) and V (</>) given by 

H'(r)=tV (0) 



(42) 



(43) 



with £ = constant. Upon considering this relation, with (38) and (39) they become equivalent and therefore we 
consider the equation (38) alone. Using above with some manipulation we obtain 



where 



Further, equations (43) implies 



and 



a 2 + U (a) = 



U (a) = f(a)- 



1 



.4 



a 2 + 1 Vor 



V{(/>) = Voe^ , (V = cons.) 



(44) 



(45) 



(46) 



(47) 



One can easily show that the following boundary condition for the dilaton is satisfied (modulo Eq. (38) and with 
reference to [1] (Eq. (38))) automatically due to the Israel junction conditions 



dcj) 2 1 dV (0) 



(48) 



At this stage we consider our solution in two different categories: a 2 ^ 1 and a 2 — 1. The latter is known also as a 
linear dilaton. 



A. The solution for a ^ 1 



In the sequel we consider the wall to be a classical, one-dimensional particle moving with zero total energy in the 
effective potential U (a) . For this purpose we employ / (a) and H (a) from the solutions (17) and (18) and to see the 
general behavior of the potential U (a) we rewrite it explicitly as 



(l + a 2 ) 2 a 2 



2 ( 2 + Q 2 ) 
P 2 Ai(^L) H« 2 
A 4 (l+a 2 ) 



U(a) 



(49) 



a 2 (l+a 2 ) 



3 + n 
l + o 



4P 2 Ai (r a \ 2 



4 / P^A 2 
3 T 



V2)-4U 1 (^)ln(^)-(32M + U 2 )(^) 



a 2 ^3 
a 2 = 3 



At this point we introduce a new parameter 7 = } + °2 (—1 < 7 < 1) and proceed to analyze the forgoing potential. 
We set also ro = 1 for convenience so that the metric function and the potential for a 2 ^ 3 (7 7^ — |) take the form 



/(a) = cj!a 1+7 + cj 2 a 1 7 + w 3 a 7 + w 4 a 1 7 
U(a) = bj\a 1+ ~< + <jJ2CL l ~ 1 + W3a -7 + u^a -1 ' 



in which 



-2Fi 



(l + 2 7 ) (I+7) 



,w 2 



(A 2 - 2P 2 A 2 ) 



8M 



,w 3 



(1 + 7) 2 



,w 4 



2P 2 Aj 
A 4 (1 + 7) ■ 



(50) 



(51) 
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while 



10 1 = uii — 



4(l + 7 ) 



2 • 



(52) 



Having the new forms of / and U, one observes that in the domain of a, for U < /, on general grounds there are 
many possibilities to be considered. We prefer to consider some cases which are not general but interesting enough 
analytically. 

The first case is to set u>i = uj 2 = for which the potential becomes simply 



U(a) = o-t (w 3 + ^) 



(53) 



Depending on the sign of Ai, M and V\ different cases may occur which are shown in Fig. 1. This figure also reveals 
that the only configuration, in this setting which admit black hole, with bouncing point are Fig. IB and IF. In Fig. 
2 we show these cases in a closer form. The bouncing points are visible and both are at 



AiP 2 (l+7) 
AM A 4 ■ 



(54) 



The other interesting setting, at which U(a) may admit two bouncing points is due to <Di = while oo 2 7^ 0. In this 
case the potential becomes 



1 



U (°) = ~1T^ (^a 2 + w 3 a + cj 4 ) 



(55) 



which clearly depends on the values of the parameters and are shown in Fig. 3. In this figure we have introduced 



A 4 7 



(c) 



46A S M 2 



(1+lfQ 2 



and 



(1 + 27) v 2 

1+7 



(56) 



(57) 



Let us add that upon setting the latter equality, Vq has no contribution in the general form of U(a) but still it changes 
the form of metric function. This is the reason that Fig.s 3B and 3C (also 3D, 3E and 3F) differ, although the settings 
are the same. In fact in these figures Vq is not the same. 



B. the Linear dilaton case a = 1 (7 = 0) 

In the case of linear dilaton with a = 1 (7 = 0), the form of the solution is given by 

~P 2 M fro 



f(r) = 2r 2 



A 4 



(7) + Hk 



AM 



together with 



H 2 



A - 



-V2 
> = 

ro J 2 \r J ' v r J \ r / \r 



(58) 



(59) 



Herein we used the condition given in Eq. (19) which in the case of k = 1 and a = 1 states 2P 2 A2 = A 2 . The DW's 
potential, therefore, reads 



U(a) = 



1 



-L-2(V 2 + V 1 )-\V 2 ) a 2 -8Ma+^ 



where r is set to one. Here also analytically one can observe that if 



P 2 \ x 



8 A 4 \A 2 



- "2d , : \ ■-\V 2 )<M 2 <^±(±,-2{V 2 + V 1 ) 



8 A 4 \A 2 



(60) 



(61) 
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and 

0<-^-2(V 2 + V 1 )-^V Q 2 , 0<Ai, (62) 
then there would be two (i.e. both maximum and minimum) bouncing points. Also if 

w(i- 2 (^ + ^)-^o 2 )>M 2 (63) 

with (62) there would be no dynamical DW universe possible (see Fig. 4) since U(a) > 0. Of course there would be 
all other case also possible which one can easily find from the closed form of the potential U(a). 

With the exception of the cases we investigated closely here many other cases behave similarly but since they should 
be treated numerically we do not study them here. For example the case of a 2 — 3 (7 = — h) can be studied only 
numerically which we ignored here. 



IV. CONCLUSION 



For a pure magnetic field in d = 4 bulk spacetime we obtain a large class of solutions in EMD theory. This class 
generalizes all previously known pure magnetic type solutions which are obtained in particular limits. The element 
that brought new extensions is to take double Liouville type coupling with dilaton in the action. The junction 
conditions induce potential on our DW which is adopted to be a 3— dimensional FRW universe. An investigation and 
plot of the intricate, induced potential reveals the possible existence of a second (maximum) bounce in our model. 
This implies an oscillatory universe on the DW between two (minimum and maximum) limits. In order to go beyond 
these bounds and give a big crunch / infinite expansion the DW universe must naturally undergo quantum tunnelling 
processes. We recall that in the Einstein - Gauss - Bonnet bulk in 5— dimensions, DWs didn't have a second bounce 
111,]. A possible extension of our model presented in this paper may be to consider electromagnetic action, study the 
Maxwell equations and find the induced charge on the DW [12j. Let us add that the absence of any DW universe 
(i.e. for the potential U(a) > 0) at all is another extreme possibility. Finally, it would be much desirable to see the 
present detailed analysis of this paper extended to higher dimensions. The difficulty originates from the fact that the 
general solution (17-18) obtained in d = 4 doesn't extrapolate to d > 4 easily. 

Acknowledgement: We are indebted much to the anonymous referee for constructive criticism and pointing out 
numerous errors in the original version of the paper. 
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Figure Caption: 

Fig. 1: Diverse plots for different sets of parameters Ai,M and 7. These range from corresponding black hole 
states (1A, IB, 1C, IE, IF) to non-black hole states (ID, 1G, 1H) in the bulk. The cases for U(a) > (ID, 1H) do 
not allow formation of DW universes at all. The figures (IB, IF) show explicit upper bounces whereas in (1A, ID, IE 
and 1G) the DW universe has infinite extension. It should be also added that cases such as f(a) < (i.e. IE outside 
the horizon and 1G every where) do not correspond to stationary spacetimes rather relate to cosmology. 

Fig. 2: With specific parameters we magnify two cases that admit upper bounces induced from the black hole 
solutions. Such a bounce enforces the universe to contract anew. Both 2 A and 2B have similar behaviour although 
the dilaton parameter differs much. 

Fig. 3: A variety of possibilities parametrized by M, 6 and Ai versus A^ c \ (For definitions, see Eq. (56) in the 
text). Depending on the sign of these parameters we display twelve different cases revealing the available richness of 
our DW universe. The bouncing from below (3B), and above (3A), or both (3D) are explicitly shown. Figures 3C, 
3J and 3L are clearly non stationary. 

Fig. 4: Various plots for linear dilaton case 7 = (or a 2 = 1). Fig.s 4A and 4B are both for non-black hole cases 
and since U(a) > they don't yield dynamic DW universes. Fig. 4C has the double bounces but doesn't correspond 
to a black hole. Fig. 4D corresponds to an extremal black hole which admits an upper bounce. 
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